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Berciu and Bhatt reply: The point raised by the 
authors of the Comment |j]] regarding the hopping term 
t(r) we use Q to describe the impurity band (IB) of 
holes in Gai_ x Mn æ As, needs clarification. A proper de- 
scription of the hopping between impurity States in such 
alloys is a very complicated and, to our knowledge, un- 
settled issue. In our work, we used a hopping parameter 
of magnitude corresponding to two isolated s-wave irn- 
purities, with t(r) < 0 in the hole representation. This 
parametrization captures two important length scales - 
the inter Mn-Mn distance, and the impurity Bohr ra- 
dius - present in a complete model of the experimental 
system. With our parameters, in the absence of mag- 
netization, we obtain an impurity band for holes whose 
density of States (DOS) is plotted in Fig.jj] (left), for a 
Mn concentration x = 0.93% using the impurity Ryd- 
berg (Ry) as a unit (lRy=112 meV for Mn in GaAs). 
The top of the hole impurity band is 3 Ry above the 
valence band, while the Fermi energy Ep for holes lies 
2.5 Ry above the top of the valence band. These appear 
to be reasonable numbers: for a random distribution of 
acceptor centers A, the top of the hole band would be 
expected to correspond to the ionization potential of the 
most stable configuration, which is likely an A 4 complex 
accomodating 4 holes in the ground impurity state. In a 
hydrogenic model, this is w 4.5 Ry ||; for centers with 
strong central cell like Mn in GaAs, it is likely reduced 
somewhat. In the mean-field approach, this is the situa- 
tion above T c , while for T < T Cl the coupling to the Mn 
spins causes the spin up and down bands to split, leaving 
the system fully polarized at T = 0 (see inset). 

Additionally, with the parameters used the IB shows a 
mobility edge close to Ep, as demonstrated by comput- 
ing the inverse participation ratio IPR=JA |<?!>(*)| 4 , where 
(j)[i) is the wavefunction amplitude at site i. For an ex- 
tended wave-function, IPR^ 1 /N r j. where TV d is the num- 
ber of sites of a finite sample. For localized wavefunctions, 
the IPR is independent of system size. The average IPR 
for the wavefunctions of our hopping Hamiltonian are 
shown in Fig. [l] (right) for Nd = 125 and 1000. The 
occupied hole States at the top of the hole impurity band 
are localized, with relatively high and size-independent 
IPR, whereas States of lower energy, which are occupied 
only at higher hole fillings, are extended with IPR’s de- 
pending on size in the expected manner. Thus, the model 
captures the proximity of the Metal-Insulator Transition, 
seen in experiment. (For x = 5%, where the magnetiza- 
tion is rnuch less anomalous, the IB top is again split a 
few Ry above the valence band, and is somewhat wider). 

On the other hånd, as pointed out in the Comment 
|lj, the choice t(r) > 0 in the hole picture, which inverts 
the IB (E — > —E), is not suitable because it leads to an 
unphysical long tail of the DOS in the gap at the densities 
of interest, and a consequent low DOS at Ep. Moreover, 
all the relevant States are extended (there is no mobility 
edge). A more realistic calculation of the IB, including 



FIG. 1. Left: DOS for hopping Hamiltonian with t(r) < 0. 
The inset shows the effect of the AFM coupling (J = 15 
meV). The er =f and er =j bands (thick lines) are 110 longer 
degenerate at low temperatures. For comparison, we plot the 
J = 0 DOS as well (thin line). Right: dependence of the 
average IPR on the energy, In both pictures x=0.0093. 

the Coulomb potential from the other Mn impurities as 
well the charged impurities responsible for the large com- 
pensation, and using the more complicated structure of 
a hole impurity wavefunction, will likely remove the un- 
physical tail and yield a DOS similar to the one used in 
our calculation, with proximity to a mobility edge. Al- 
ready the tail is absent in e.g. Ref. 3, for a hydrogenic 
lattice. Since the nature of compensation in GaMnAs is 
still an open question, we opted to use a simple model 
that gives a physically acceptable description of the IB. 

Our goal in Ref. [|| was to use this simple IB model 
to point out non-trivial effeets of disorder in Mn po- 
sitions on the shape of the magnetization curve M(T ) 
and the critical temperature Tq. To our knowledge, all 
previously published studies had neglected this aspect. 
Studies appearing since support our claim of inereased 
Tc with inereased Mn disorder Q. While we agree that 
a better modeling of the impurity band and inclusion 
of the valence band and other factors such as screen- 
ing (on a proper, local scale, taking into account strong 
charge inhomogeneities) are necessary to achieve a proper 
quantitative description especially in the metallic regime 
x > 0.03, we maintain that the underlying physics cap- 
tured by our simple model of Ref. [|| is essentially correct. 
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